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Q\ ■ Abstract 

*-£h \ Mixed dark matter scenario can reconcile the COBE data and the observed large 

i" 1 ■ scale structure. So far the massive neutrino with a mass of a few eV has been the 

^h. only discussed candidate for the hot dark matter component. We point out that 

c ! ! the hadronic axion in the so-called hadronic axion window, f a ~ 10 6 GeV, is a 

perfect candidate as hot dark matter within the mixed dark matter scenario. The 
current limits on the hadronic axion are summarized. The most promising methods 
to verify the hadronic axion in this window are the resonant absorption of almost- 
monochromatic solar axions from Ml transition of the thermally excited 57 Fe in the 
Sun, and the observation of the "axion burst" in water Cerenkov detectors from 
another supernova. 
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1 Introduction 



The cold dark matter (CDM) dominated universe with scale-invariant primordial density 
fluctuation has been the standard theory of structure formation. After COBE has found 
the finite density fluctuation in the cosmic microwave background radiation (CMBR), 
the standard CDM scenario was found to give too much power on smaller scales. Many 
modifications to the standard CDM scenario were proposed which solve the discrepancy: 
by introducing a small Hot Dark Matter (HDM) component [|TJ, by "tilting" the primor- 
dial density fluctuation spectrum 0, by assuming a finite cosmological constant B, or 
by introducing particles (such as v T ) whose decay changes the time of radiation-matter 
equality 0. At this point, there is no clear winner among these possibilities.^ 

In this letter, we revisit the mixed dark matter (MDM) scenario from the particle 
physics point of view. This scenario has attracted strong interests because there has 
been a natural candidate for the HDM component: massive neutrino(s). A neutrino 
with a mass of a few eV can naturally contribute to a significant fraction of the current 
universe. However, it has not been easy to incorporate the HDM together with other 
neutrino "anomalies," unless all three generation neutrinos (possibly together with a 
sterile neutrino) are almost degenerate, and their small mass splittings explain various 
"anomalies." Such a scenario may be viewed as fine-tuned. Especially, the atmospheric 
neutrino anomaly is quite significant statistically now thanks to the SuperKamiokande 
experiment, which suggests the mass squared difference of Am 2 = 1CT 3 — 1CT 2 eV 2 between 
the muon and tau neutrinos. If we view the situation from the familiar hierarchical fermion 
mass matrices, it suggests the tau neutrino mass of 0.03 - 0.1 eV, and it appears difficult 
to accommodate the HDM based on massive neutrinos. 

We point out that the hadronic axion |7J can be an alternative motivated candidate 
for the HDM component in the MDM model. Axion has been proposed as a solution 
to the strong CP problem in the QCD, and the hadronic axion (or KSVZ axion) is one 
version which predicts small coupling of the axion to the electron. There has been known 
a window of f a ~ 10 6 GeV allowed by existent astrophysical and cosmological constraints 
if the axion coupling to photons is suppressed accidentally. This is referred to as the 
"hadronic axion window." Our main observation is that this window gives exactly the 
right mass of m a ~ a few eV and the number density of the axion appropriate for the 
HDM component in the MDM scenario. 

2 Hadronic Axion 

First, let us review the hadronic axion model J7[. The most important feature of the 
hadronic axion is that it does not have tree-level couplings to the ordinary quarks (u, d, 

1 However, a large "tilt" is difficult to obtain in many inflationary models. tCDM can be tested well 
by £?-factory experiments in the near future The recent data from high-redshift supernovae prefer 
ACDM H , but the possible evolution of supernovae needs to be excluded by more systematic comparison 
between nearby and high-z supernovae. 



1 



s, c, b, t) and leptons (e, u e , fi, r, v T ). In this framework, we introduce new fermions 
which have Peccei-Quinn (PQ) charges, while ordinary fermions do not transform under 
U(1)pq. Some of those new fermions, which we call PQ fermions hereafter, also have 
SU(3)c quantum numbers. After the PQ symmetry is broken spontaneously, axion a 
appears as a pseudo-Nambu-Goldstone boson of the PQ symmetry. 
The axion a couples to the photon with the operator 

£a77 T^Qa-yydt^ 1 ^ FuvFpa = Z ^ FuvFp a , (1) 

8 167T f a 

where f a is the axion decay constant. This interaction is induced by the mixing to the 
light mesons (tt°, rj, rj', and so on) as well as by the triangle anomaly of the PQ fermions. 
By using the chiral Lagrangian based on flavor SU(2) L x SU(2)#, we can estimate the 
coefficient C a77 as || 



E PQ 2(4 + z) 
a ~ n N 3(1 + 2)' 



where z = m u /md which is estimated to be 0.56 by the leading order perturbation in 
quark masses in the chiral Lagrangian. (Hereafter, we use z = 0.56 for our estimation, 
unless we discuss quantities which are sensitive to the uncertainty in z.) In Eq. (0), the 
first term is from the U(l) em anomaly of the PQ fermions, while the second term is due to 
the mixing between axion and light mesons. Simultaneously, we also obtain the formula 
for the axion mass as 

m a = ^ ~ 6.2 eV x (/ a /10 6 GeV)-\ (3) 

± r Z J a 

where f n ~ 93 MeV is the pion decay constant, and is the pion mass. 

With this axion-photon-photon coupling, axion decays into two photons with the life- 
time 



a 2 C 2 _3 



T77' 



2567T3 fl 



1 

~ 1.2 x 10 12 yr x C~± K/10 eV)~ 5 . (4) 



Notice that the lifetime of the axion is longer than the age of the Universe for m a ~ 10 eV 
and C a77 ^ 1, and hence primordial axions are still in the Universe. However, as we will 
see later, radiative decay of the axion may affect the background UV photons in spite of 
the long lifetime. 

Here, we comment that C a77 is significantly affected by uncertainties in the chiral 
Lagrangian with which the mixing effect is usually calculated. First of all, the accuracy 
of the SU(2)^x SU(2)h chiral Lagrangian is tested up to about 5 - 10 %. For example, 
by using the pion decay constant estimated from the leptonic decay width of 7r ± , r(7r° — > 
7 + 7) is calculated to be 7.73 eV 0, while experimentally, it is measured to be 7.7 ± 
0.6 eV [0. (Even though the center value given in Ref. [0] is in a good agreement, 
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the single best measurement suggests the width to be 7.25 ± 0.23 eV |Tl|, which is about 
6 % off from the chiral Lagrangian prediction.) Furthermore, /„. estimated from the 
process e + + e~ — > 7r° + e + + e~ |12| is about 10 % smaller than the one from the 
leptonic decay of ir^ 1 ||. Therefore, we may expect 5 - 10 % error in the calculation 
of the mixing effect from chiral Lagrangian. Another uncertainty is from the so-called 



Kaplan-Manohar ambiguity |T3[ . Within the lowest-order chiral perturbation theory, z is 
estimated to be 0.56. However, under the SU(2)lX SU(2)# flavor symmetry, the quark 
mass matrix M = diag(m tt , m^) and its conjugate (ia 2 )M*(— ia 2 ) = diag(m^,m*) have 
the same transformation properties, and hence the following shifts are allowed: m u — > 
m' u = m u + em*,, — > m' d = m^ + em*, where e is an unknown parameter |T3. Since the 



parameter e is arbitrary, z = m u jmd cannot be determined from the meson masses alone.Q 
In particular, z much smaller than 0.56 (or even z = 0) may be allowed if we take this 



ambiguity into account [13[]. This ambiguity cannot be resolved based on meson masses 
only, but can be by using the baryon masses to some extent. The uncertainty, however, 
remains large The mixing contribution to C ail is affected by this uncertainty in z. 

As we will see later, C ail is constrained to be less than 0.01 - 0.1 from astrophysical 
arguments for the axion decay constant we are interested in. In general, C ail 1 is 
possible if we adopt an accidental cancellation. With the lowest order chiral Lagrangian, 
cancellation occurs when Epq/N = 2(4 + z)/3(l + z) ~ 1.95, but this estimation may not 
be so reliable. We believe that a better understanding of the quark masses is necessary 
to pin down the value of Epq/N for the accidental cancellation. With the current best 
knowledge, it is clear that the cancellation is quite possible for models with E-pq/N ~ 2 
if we take the effects we discussed above into account. In particular, the possibility of the 
value obtained in grand- unified theories (Epq/N = 8/3) may not be excluded. 

The axion is also coupled to fermions: C a ff = gaffdfi^f, which can again be esti- 
mated by using the chiral Lagrangian. Importantly, the hadronic axion does not couple to 
ordinary quarks and leptons at the tree level. Therefore, in particular, the axion-electron- 
electron coupling has an extra loop suppression factor ||15|| : 



3a 2 m e (E PQ 2(4+^) 1 

9aee = ^ j~ < — Hf a /m e ) - ^ - ^ ln(A QCD /m e )| . (5) 

On the other hand, mixing effects induce an axion-nucleon-nucleon coupling, even though 
the axion-quark-quark coupling vanishes at the tree level for a hadronic axion: 

9aNN = ^\(Fao T Faz);^— + (F A0 ± F A3 ) Z ) , (6) 
fa { 2(1 + z) 2(1 + z) J 

where ~ 940 MeV is the nucleon mass, and upper (lower) sign is for neutron (proton). 
The axial- vector isovector contribution has been quite well understood to be ~ —1.25 

2 In the SU(3)lxSU(3)# chiral Lagrangian, the effect is formally higher order in quark masses, and 
hence e ~ m s /(47r/ ff ). Still, the ambiguity in z is rather large. 

3 However, if z = 0, strong CP problem is solved without introducing an axion. Therefore, we do not 
consider this possibility in this letter. 
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from the neutron /3-decay. Isoscalar part Fao used to be more ambiguous, since this 
quantity depends on the flavor-singlet axial-vector matrix element S (with S = Au + 
Ad + As in Ref. [fL7|0 as Fao — — 0.675 — 0.20, where the constant piece is determined 



by the hyperon /3-decay. In Ref. [|lYj], however, S was estimated from experimental data 
including higher order QCD corrections, resulting in S — 0.27 ± 0.04. Even though 
possible systematic uncertainties are not included in this calculation, we use this result as 
a reference when we discuss axion-nucleon-nucleon coupling. Because of these interactions, 
f a is constrained by the axion emission from SN1987A. 



3 Constraints on Hadronic Axion 

Next, we summarize the constraints on the hadronic axion. In the later discussion, we 
will be interested in the case of f a ~ 10 6 GeV so that hadronic axion becomes a good 
candidate of the HDM. Therefore, in this section, we pay an attention to this case. 

Most importantly, the coupling of the hadronic axion to the electron is loop suppressed, 
as can be seen in Eq. (H). Therefore, the constraint on the axion-electron-electron coupling 



from the cooling of the red giant []18 , [19| can be evaded. One can compare the current 
best upper limit (g aee ^ 2.5 x 10~ 13 |T^]) with Eq. (|5p, and see that g aee for f a ~ 10 6 GeV 
is smaller than the bound from the red giant for values of E-p^/N < 7. 

A non-trivial constraint comes from the emission of the axion from a supernova. If 
an axion couples to nucleons strongly, the axion can be produced in the core of the su- 
pernova, and the axion emission may affect the cooling process of the supernova. In 
particular, the Kamiokande group and the 1MB group measured the flux and duration 
time of the neutrino burst emitted from the SN1987A, and their results are consistent 
with the generally accepted theory of the core collapse. Therefore, they confirmed the 
idea that most of the energy released in the cooling process is carried off by neutrinos. 
If axion carries away too much energy from the supernova, it would conflict with those 
observations. The axion flux from the supernova can be suppressed enough in two pa- 
rameter regions. If axion-nucleon-nucleon interaction is weak enough, the axion cannot 
be effectively produced in the core of the supernova. Quantitatively, for f a > 10 9 GeV, 
the axion flux can be small enough not to affect the cooling process |2(J . On the contrary, 
if the axion interacts strongly enough, the mean free path of the axion becomes much 
shorter than the size of the core, and hence the axions cannot escape from the supernova. 
In this case, axion is trapped inside the so-called "axion sphere," and the axion emission 
is also suppressed. (In this case, axions are emitted only from the surface of the axion 
sphere; this type of the axion emission is often called "axion burst.") Quantitatively, for 
f a < 2 x 10 6 GeV (or equivalently, m a ^ 3 eV), the axion luminosity from SN1987A is 
suppressed enough 

For f a ^ 2 x 10 6 GeV suggested from the cooling of supernova, we have another 
constraint from the detection of axions in water Cerenkov detectors. In this parameter 
region, axion flux from the axion burst is quite sizable for its detection, even though it does 
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not affect the cooling of SN1987A. If the axion-nucleon-nucleon coupling is strong enough, 
axions may excite the oxygen nuclei in the water Cerenkov detectors ( 16 + a — > 16 0*), 
followed by radiative decay(s) of the excited state. If this process had happened, the 
Kamiokande detector should have observed the photon(s) emitted from the decay of 16 0*. 
Due to the non-observation of this signal, / a <3x 10 5 GeV is excluded [f21~1 . 



Another class of constraint is from the axion-photon-photon coupling. Because of 
this coupling, axion can be produced in Primakoff process in the presence of external 
electromagnetic field, and it also decays into two photons, which result in constraints on 
the (model-dependent) axion-photon-photon coupling. 

One of the important constraints comes from the cooling of the horizontal-branch (HB) 
stars. If the axion-photon-photon coupling is too strong, axions are produced in the HB 
stars through the Primakoff process, and the emission of the axions affects the cooling 
of the HB stars. Then, the lifetime of the HB stars becomes shorter than the standard 
prediction, and the number of the HB stars are suppressed. However, the number of the 
HB stars are in a good agreement with theoretical expectations, and hence we obtain the 



upper bound on the axion-photon-photon coupling | 22| : 



g ail < 6 x l(T n GeV . (7) 

The important point is that g ail has two sources: the electroweak anomaly of the PQ 
fermions and the mixing between the axion and light mesons (see Eqs. ([I]) and (0)). 
Furthermore, the mixing effect is usually calculated by using the chiral Lagrangian, and 
there is some uncertainty as discussed earlier. Therefore, it is difficult to convert the 
constraint (|7|) to the constraint on the PQ scale f a . In fact, due to the model dependence, 
we only have an upper bound on the coefficient C a77 : 

C a77 < 0.05 x (/ a /10 6 GeV). (8) 

Notice that, in principle, any value of f a can be viable with the cooling of the HB stars, 
if we adopt an accidental cancelation in C a77 . 

Another important constraint is from the effects of the radiative decay of the axion 
on the background UV photons. As noted in Eq. (Jl|), axion is coupled to photons, and 
it decays into two photons with the lifetime given in Eq. (f|). Even though the lifetime 
is longer than the age of the Universe, some fraction of the axion decays and we may see 
the emission line. 

Constraint from the UV extragalactic light is discussed in Refs. [jR], ^3|, Since 
the lifetime of the axion is longer than the age of the Universe, intensity of the photon 
is proportional to the inverse of the lifetime. Therefore, the intensity becomes smaller as 
the axion-photon-photon coupling gets weaker, and non-observation of the signal sets an 
upper bound on C a77 . Overduin and Wesson looked for the emitted photon from the axion 
in the extra galactic light, and no signal of the axion was found. From their observation, 
they derived the upper bound on C a77 of 0.72 (m a = 3.8 eV) to 0.014 (m a = 13.0 eV) ||. 

More stringent constraint may be obtained if we observe the photons emitted from 
the axions in clusters of galaxies. At the center of a cluster, axions are expected to be 
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Figure 1: Astrophysical constraints on the hadronic axion model from the cooling of the 
supernova, axion burst, cooling of the HB stars, extragalactic light [^2] (square), and 
emission line in clusters of galaxies [Q (triangle). Shaded region is excluded, and C a77 
larger than squares and triangles are inconsistent with observations for fixed value of f a . 



gravitationally trapped, and its density is more enhanced than the cosmological density. 
Therefore, the emission lines may be more intense than the one from the extra galactic 
sources, and the constraint may be more stringent. With three samples of clusters, Ressell 
obtained the upper bound on C ail of 0.12 (m a = 3.5 eV) to 0.008 (m a = 7.5 eV) ^3|, which 
is about one order of magnitude more stringent than the constraint from extra galactic 
background light. However, it is possible that the lines of sight of the particular galactic 
clusters are obscured by absorbing material, resulting in too stringent constraint p4j. If 
we adopt this argument, this constraint may be evaded. 

All the constraints mentioned above are summarized in Fig. [I]. As we have discussed, 
the hadronic axion with the axion decay constant in the following range is still viable with 
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all the astrophysical constraints (if C a77 is small enough): 

3 x 10 5 GeV < f a < 2 x 10 6 GeV (20 eV > m a > 3 eV). (9) 

Notice that the constraints based on the axion emission from SN1987A is relatively model- 
independent. That is, in the hadronic axion model, the axion-nucleon-nucleon coupling 
is from the mixing between the axion and the light mesons, and hence it is independent 
of the U(1)pq charges of the PQ fermionsfj 

Finally, we comment on the constraint from the cooling of the red giants and the 



HB stars due to the axion-nucleon-nucleon coupling [2S|. The axion-nucleon-nucleon 
coupling would allow an axion emission from red giants and the HB stars, and cause an 
additional energy loss rate which is proportional to m? a . This extra energy loss changes 
the brightness of these stars, and it also modifies the relative numbers of the red giants 
to the HB stars. Observed values of these quantities are in reasonable agreements with 
theoretical calculations, and hence we can obtain the upper bound on the axion emission 
rate. The constraint is quite sensitive to the flavor- singlet axial-vector matrix element 
S = Au + Ad + As, since axion-nucleon-nucleon coupling depends on S. For S = 0.27 as 
suggested in Ref. [fF7|, axion mass smaller than about 12 eV is still allowed,^ and larger 



axion mass is still viable if we adopt sizable uncertainty in S [J25|. Therefore, we concluded 



that most of the parameter region for the axionic HDM is still alive. 



4 Thermal Relic of Hadronic Axion 

We have seen in the previous section that the hadronic axion with the decay constant in 
the window 3 x 10 5 GeV ^ /„ ^ 2 x 10 6 GeV is astrophysically allowed as long as the 
axion-photon-photon coupling is sufficiently small. Now, we are in the position to discuss 
how the hadronic axion can be a good candidate for HDM. For this purpose, remember 
that the relevant mass range for the HDM is 1 eV - 10 eV, corresponding to the PQ scale 
of f a ~ 10 6 GeV (see Eq. (^)), if the axion decouples around the same stage as when the 
neutrinos do. 

For f a ~ 10 6 GeV, the most important source of the primordial axions is the thermal 



production, rather than the coherent oscillation [16, [L5[. Because of the couplings to 



nucleons (and to pions), axion are thermalized when T > 30 — 50 MeV for f a ~ 10 6 GeV. 



In the most recent calculation |15[], the axion density is estimated as [p a / Pv]t~i MeV 
0.4 — 0.5, with p a (p u ) being the energy density of the axion (neutrino of one species), or 
equivalently, 

Tt 

— ~ 0.02, (10) 
s 

4 It does suffer from the uncertainty in z mentioned earlier, however [ p"5| . 

5 The authors of Ref. |2j| used F A0 = -0.675 - 0.23 from hyperon (3 decay without SU(3) breaking 
effects. A direct measurement, however, suggests — 0.675 — 0.20 |2(J, and makes the S in their plot 
effectively smaller by 0.04. 
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where n a is the number density of axion, and s is the total entropy density. (Here, we 
used [pa/Pulr^i Mev — 0.45.) Then, the relic density of the axion is given by 

771 71 

n a = ~ 0.2 x h£(m a /10 eV), (11) 

Pc 

where /150 is the Hubble constant in units of 50 km/sec/Mpc. Thus, for m a ~ 10 eV, Q a 
can be 0.1 - 0.2 which is the requirement for the HDM in the MDM scenario. Note that 
the hadronic axion discussed here is a thermal relic with its mass of ~ 10 eV. Therefore, 
the axion here is a relativistic particle when the galactic scale crossed the horizon, and 
behaves as HDM.0 

Comparing with Eq. (pTTj) , the window (|9]) is exactly where the axion has the right 
mass and number density to be the HDM component in the MDM scenario. 

One may worry about the effect of the hadronic axion on the big-bang nucleosynthesis 
(BBN). At the time of the BBN, energy density of the axion is sizable ([p a / Pv]t~i Mev — 
0.4 — 0.5), and it raises the freeze out temperature of the neutron by speeding up the 
expansion rate of the Universe. As a result, in our case, more 4 He is synthesized than in 
the standard BBN case flT5|. A few years ago, the observed value of the primordial 4 He 



abundance seemed to be unacceptably smaller than the theoretical prediction ||29|| . If this 



was true, a hadronic axion with f a ~ 10 GeV could be extremely disfavored. However, 
the current situation is more controversial. Recently, both for D and 4 He, several new 
measurements have been done to determine their primordial abundances, but the results 
are not consistent with each other; some group reports low D abundance while the other 
results are much higher, and the same for 4 He. In particular, if we adopt a high value 
of the observed 4 He abundance our scenario is consistent with the BBN. Since it 
is too premature to judge which measurements are reliable, we do not expect any solid 
argument based on the BBN which rules out the hadronic axion as the HDM component 
in the MDM scenario. 



5 Prospect for Detecting Hadronic Axion 

So far, we have seen that the hadronic axion in the current allowed parameter range almost 
automatically becomes appropriate for HDM. As discussed, this scenario is consistent 
with all the astrophysical constraints, if the axion-photon-photon coupling is suppressed 
enough, presumably by an accidental cancellation. 

6 It is interesting to note that the axion decay constant required in this scenario is rather close to the 
so-called messenger scale in models with gauge mediation of supersymmetry breaking |27j , as well as the 
mass scale of the right-handed neutrino in the sneutrino CDM scenario [Q . It is conceivable that the 
field S which generates the supersymmetric and supersymmetry-breaking masses of the messengers carry 
the PQ charge and the messengers are the PQ fermions. The same field S can generate the required size 
of the right-handed neutrino mass in the sneutrino CDM scenario. The original scale of supersymmetry 
breaking, however, needs to be raised to make the gravitino heavier than the sneutrino, which can be 
achieved by making the messenger U(l) coupling constant somewhat small, ~ 0.03. 
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However, this scenario can be tested in the future in several observations. One possibil- 
ity is to use the observation of the diffuse background UV photon. Accuracy of the current 
observation just excluded the axion-photon-photon coupling down to C ary ^ 0.1 — 0.01, 
as we have discussed. However, if the background photon spectrum will be well measured 
with a better resolution, the emission line from the axion decay may be found in the 
background photon spectrum. However, as we emphasized, C ail is a model-dependent 
parameter. Therefore, a non-observation of the signal cannot exclude the possibility of 
hadronic axion HDM definitively, because of a possible accidental cancellation in C ail . 

Therefore, a detection of hadronic axion which does not rely on axion-photon-photon 
coupling is strongly favored. One such possibility is to detect an axion burst from a future 
supernova at SuperKamiokande (or, in general, water Cerenkov detectors). An important 
point is that newer water Cerenkov detectors (like SuperKamiokande) have much larger 
fiducial volume than Kamiokande, and hence we can expect a larger event rate. Therefore, 
a hadronic axion with f a ~ 10 6 GeV can be tested with a future supernova of the size 
and the distance of SN1987A, even though SN1987A could not exclude this possibility. 

Calculation of the event rate suffers from the uncertainties in the axion-nucleon scatter- 
ing cross section and modeling of supernovae. However, the detection of the signal appears 



plausible. For example, by rescaling the result given in Ref. [21|, we expect a few events 
at SuperKamiokande for a supernova of the same size as SN1987A for f a ~ 10 6 GeV. Of 
course, if a new supernova will be closer than SN1987A, we can expect larger number of 
events, and the hadronic axion HDM can be tested much easier. 

Another interesting novel idea is due to Moriyama |31| . In the Sun, thermally excited 
57 Fe nuclei can decay by emitting axions. Thanks to the Doppler broadening of the axion 
energy due to the thermal motion of 57 Fe, the same nuclide can resonantly absorb the 
axion. The detection rate was estimated and can be as high as 1 event/day/kg or more. 
A search was already performed along this line |E2J even though they used a small target 
of 0.03 g to detect 14.4 keV gamma-ray escaping the target rather than the bolometric 
method suggested. They obtained an upper bound on the axion mass of 745 eV. Another 
experimental effort to detect solar axions is underway and may reach the axion mass as 
small as 3 eV in a few years [B3 . 



6 Conclusions 

In this letter, we have pointed out that the hadronic axion in the hadronic axion window 
(f a ~ 10 6 GeV) can automatically be a good candidate of the Hot Dark Matter component 
in the mixed dark matter scenario. In order to evade an astrophysical constraint from 
the background UV light, axion-photon-photon coupling has to be suppressed in the 
hadronic axion window, probably by an accidental cancellation. This scenario may be 
tested by detecting the axion burst from a future supernova in water Cerenkov detectors, 
or detecting solar axions using resonant absorption. 



9 



Acknowledgements 



We thank Shigetaka Moriyama for letting us know his work. HM thanks Gia Dvali, Markus 
Luty, Joe Silk, and Mahiko Suzuki for useful conversations. This work was supported in 
part by the U.S. Department of Energy under Contracts DE-AC03-76SF00098, in part by 
the National Science Foundation under grant PHY-95-14797. HM was also supported by 
Alfred P. Sloan Foundation. 



References 

See, for example, J.R. Primack, |astro-ph / 9 70 7285) . 

R. Davis et al, Phys. Rev. Lett. 69, 1856 (1992); A. Liddle and D. Lyth, Phys. Lett. 
B291, 391 (1992); F. Lucchin et al, Astrophys. J. 401, L49 (1992); J. Gelb et al, 
ibid. 403, L5 (1993); R. Cen et al, ibid. 399, Lll (1992). 

J. Ostriker and P. Steinhardt, Nature 377, 600 (1995). 

S. Dodelson, G. Gyuk, and M. S. Turner, Phys. Rev. Lett. 72, 3754 (1994); M. White, 
G. Gelmini, and J. Silk, Phys. Rev. D51, 2669 (1995); S. Bharadwaj and S. K. Kethi, 
|astro-ph/9707143 . 



J.L. Feng, T. Moroi, H. Murayama and E. Schnapka, |hep-ph /9 7094 1 1| , to appear in 
Phys. Rev. D. 

S. Perlmutter et al, Nature 391, 51 (1998). 

J.E. Kim, Phys. Rev. Lett. 43, 103 (1979); M.A. Shifman, A.I. Vainstein, and V.I. 
Zakharov, Nucl. Phys. B166, 493 (1980). 

See, for example, EW. Kolb and M.S. Turner, The Early Universe (Addison- Wesley, 
1990). 

W.J. Marciano and A. Sarlin, Phys. Rev. Lett. 71, 3629 (1993). 

Particle Data Group, R.M. Barnett et al, Phys. Rev. D54, 1 (1996). 

HW. Atherton et al, Phys. Lett. B158, 81 (1985). 

H.-J. Behrend et al, Z. Phys. C49, 401 (1991). 

D.B. Kaplan and A.V. Manohar, Phys. Rev. Lett. 56, 2004 (1986). 

M.A. Luty and R. Sundrum, Phys. Rev. D52, 5202 (1995). 

S. Chang and K. Choi, Phys. Lett. B316, 51 (1993). 



10 



[16 

[ir; 
[is; 

[19 

[20; 

[21 

[22; 

[23; 
[24 
[25; 
[26; 
[27 

[28; 
[29; 

[30 

[31 

[32; 



M.S. Turner, Phys. Rev. Lett. 59, 2489 (1987). 
J. Ellis and M. Karliner, |hep-pri/9601"280| . 

G.G. Raffelt, Phys. Rev. D33, 897 (1986); D.S.P. Dearborn, D.N. Schramm, and G. 
Steigman, Phys. Rev. Lett. 56, 26 (1986); T. Altherr, E. Petitgirard, and T. del Rio 
Gaztelurrutia, Astropart. Phys. 2, 175 (1994). 

G. Raffelt and A. Weiss, Phys. Rev. D51, 1495 (1995). 

M.S. Turner, Phys. Rev. Lett. 60, 1797 (1988); A. Burrows, M.S. Turner, and R.P. 
Brinkmann, Phys. Rev. D39, 1020 (1989); A. Burrows, M.T. Ressell, and M.S. 
Turner, Phys. Rev. D42, 3297 (1990). 

J. Engel, D. Seckel, and A.C. Hayes, Phys. Rev. Lett. 65, 960 (1990). 

G.G. Raffelt, Stars as Laboratories for Fundamental Physics: The Astrophysics of 
Neutrinos, Axions, and Other Weakly Interacting Particles (University of Chicago 
Press, 1996). 

M.T. Ressell, Phys. Rev. D44, 3001 (1991). 
J.M. Overduin and PS. Wesson, Astrophys. J. 414, 449 (1993). 
W.C. Haxton and K.Y. Lee, Phys. Rev. Lett. 66, 2557 (1991). 
S. Y. Hsueh et al, Phys. rev. D38, 2056 (1988). 

M. Dine, A. Nelson, and Y. Shirman, Phys. Rev. D51, 1362 (1995); M. Dine, A. Nel- 
son, Y. Nir, and Y. Shirman, Phys. Rev. D53, 2658 (1996). 

L.J. Hall, T. Moroi, and H. Murayama, [hep-ph/9712515| , to appear in Phys. Lett. B. 

N. Hata, R.J. Scherrer, G. Steigman, D. Thomas, T.P. Walker, S. Bludman, and 
P. Langacker, Phys. Rev. Lett. 75, 3977 (1995); C.J. Copi, D.N. Schramm, and M.S. 
Turner, Phys. Rev. Lett. 75, 3981 (1995). 

Y.I. Izotov, T.X. Thuan, and V.A. Lipovetsky, Astrophys. J. Suppl. Series 108, 1 
(1997). 

S. Moriyama, Phys. Rev. Lett. 75, 3222 (1995). 

M. Krcmar, Z. Krecak, M. Stipcevic, A. Ljubicic, D. A. Bradley, IRB-ZEF-98/01, 
|nucl-ex/980i005 . 



[33] S. Moriyama, private communications. 



11 



